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ABSTRACT 


This work 1S an application of the methods in obtaining 
the correlation functions. In particular, the method of 
correlation without a "pure" time-delay 1s presented 
together with the concept of “orthogonal filters", which are 
Laguerre function type filters. Of these filters, the 
non-symmetric Laguerre type 1s analyzed and used to realize 
a practical correiator designed for low frequency Signals. 

The correlator WaS computer-Simulated by the DSL 
subroutine and the results of the autocorrelation of a 155 
HZ sine wave were compared to the results obtained for the 
autocorrelation of a Similar wave in the actual correlator. 
A detail description of the design of the correlator and of 


the DSL program used are also presented. 
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Te INTRODUCTION 


The detection of a Signal in the presence of noise 
requires a process which becomes more complex when little is 
known of the noise characteristics and of the signal. To 
this end, correiation analysis 1s one of the most valuable 
technigues used to filter noise from very low level signals. 
Pimeaceye ene PLOplecn Of getting the signal out of noise is 
the major function of analyzing data in many areas today. 

When properly used and understood, correlation provides 
the engineer with a technique as fpowerful as Fourier 
Analysis or any of the other classical techniques. 
Essentially, correlation analysis and Fourier transform can 
be thought of as duais of each other Since in general the 
transform of any extended signal in time 1S a narrow 
frequency signal and viceversa, and the autocorrelation of a 
Single frequency sine wave extends over the entire tine 
Getay axis. Femilisrization then with this technique and 
the ability to shift from the time to frequency to delay 
axeS in order to extract more and more useful information 
from a multitude of signals is of the most importance for 
communications and related fields. 

The mathenatical background of the correlation 
functions, autocorrelation and ECEOGSSECOFECIaGt LON, is 
presented first ype et . eo | together with Chews 
properties and Fourier transform correspondence. Then, the 
method of obtaining the correlation functions are briefly 
explained. in particular, the practical application of the 
correlaticn method without a "pure" time delay is presented. 
This method is based on the application of the concept of 
orthogonal fame Ss is which are lLaguerre Functions’ type 
Emvters. 

A practical application of this method was carried out 
in the laboratory to verify the analysis. The analysis was 


done by means of the DSL subroutine which simulated the 





correlator. The experimental results of the autocorrelation 
of a sine wave of 155 HZ2 were as predicted by the 


Ssinulation. 


1G 





Ite CORRELATION ANALYSISe 


ee Se See a ee == > a eee ce 


Correlation analysis is a technique used to determine 
the spectral characteristics of a signal or the similarity. 
of two different Signals. 

A good method of measuring the similarity between two 
Signals is to multiply them together, ordinate by ordinate 
and to add the products over the duration of the waveforms. 
The result is a single number which represents the 
Similarity between the two wave forms.[Ref. 3] 

The two correlation functions, autocorrelation and 
crosscorrelation, can be explained as follows (Ref. 1,2]: 

Suppose that a stationary, physical process is producing 


meiee time functions ae) Bele xX (bt) Samultaneously. 
n 


(A process is said to be stationary if its amplitude 
statistics do not change with time). It is also assumed 
that: 

(1) The time functions being generated are not zero. Or 
they have zero mean if the time functions correspond to a 
stochastic precess 

(2) They do not have a DC component. 

(3) They can be simple or complex periodic waves. 

(4) They may vary in a random fashion. 

Under these conditions, the two correlation functions 


are defined as: 
Ae AUTOCORRELATION FUNCTION®e [kef. iy et 


The tine average autocorrelation function for 
deterministic or known Signals is defined for two different 
classes of signals: 

a) Signals of finite average power. 

b) signals of finite energy or pulse signals 


For signals of finite average power the definition 1s: 
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ik 
K u) = lin 1/2T \ X (t)x (ttuddt 2204 
fl) = die 7 , Ct) x, (eta) (2.01) 


where "u" represents a time-shift. 


For signals of finite energy the definition is: 


CO 


R (u) = \ x (t) x (ttu) dt (2.02) 
XX “ee 1 1 


The autocorrelation of functions generated by a 


stochastic process is defined as: 


R(u) = Ef a cose: 7} (23,03) 


and if the process is ergodic, the autocorrelation function 


is defined as: 


R  (u) = R(u) (2.04) 
XX 


In general, the autocorrelation function of a waveforn 


ee 1S a graph of the similarity between the waveforn 
ae) and a delayed replica of itself, Ra as a function 


of the time shift "u", 

Based on the mathematical definition OF 
autocorrelation, Figure 1 shows an elementary device 
designed to obtain the autocorrelation function. Here, the 
Signal is multiplied by its delayed replica and the result 
1S averaged over a sufficiently long time. The final result 
is the autocorrelation function if u has also been varied 
during the process. This 1s essentially a time-delayed 


autocorrelator. 


12 
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1. Properties of the Autocorrelation Function. 
The final output of figure 1, that is, the 
autocorrelation function, will always have the following 
Characteristics for any x(t): 
(1) When u = OQ the product is maximum. This is 


represented mathematically by Schwartz Inequality: 


IR (uy tl << R_ (0) (2.05) 
xXx xX 


(2) The value of the output at u = 0 represents the 


total power of the Signal. [If eee is a voltage, then for 


2 
m= 0 , R (0) = a) is the average power of the signal 
XX 


(as measured on a 1 ohm resistor). 
(3) The autocorrelation function is also a 
memetion of 


n 


(4) The Shape of the function is characteristic 
Gieene OLLginai Signal eee: 


(5) The autocorrelation function is an even 


mumerion of u. HMathenatically: 


R (u) = R- (-u) (2.06) 


(6) its nae contains periodic frequency 

components, then R (uj) will contain each of these freguency 
ek 

components. This means that for certain time functions, if 


Be) is given then R (u) 1s known which is useful in 
xX 


detecting signals py radar. 
(7) For a real process of this kind, the 
autocorrelation function approaches zero as u -->@ .. In 


other words, the signal loses coherence as the delay u is 


14 





increased. So the coherence time of the original signal is 
measured by that value of u that produces a significant 
reduction in the value of the autocorrelation function. 
The coherence time of any process producing very wideband, 
uniforn (white) aoise, is practically zero or very short 
because its instantaneous value is nearly independent of the 
value at any other time. 

(8) The autocorrelation function and the power 
spectral density are Fourier transform pairs. So they 
contain the same information, nowever, the autocorrelation 
function contains this information in the form of a function 


of time rather than frequency. 


Mathematically: 
a -jwu 
g(w) = Va (u)e du (2.07) 
Vee se 
a Jwu 
R (u) = 720 \ p(w) « dw (2210 5) 
xX — 


The power spectral density of finite average power 


Signals is: 


on 
= lim 1/2T X X Dae 
Z (w) ae f2T ana a (Ww) ( ) 
2 
= lim 1/27 |X (w)] (2210) 
Tyo 1T 


where * means the conjugate value of the quantity it refers 
to. The signal has been observed through a window of 
observation of duration fT. 

The power spectral density for signals of finite 
energy 1S: 


na Z 
Bw) = Semen SS (w) = ee (2.11) 


iS 





This is also true for stochastic processes in which 
case: 


R (u) <=> g(w) (2.12) 
XxX 


and 


2 
= lin 1 . 
B(w) aon J 20 Be ee } (2.13) 


If the Fourier components of the waveform are 
Squared in amplitude, set into phase at the origin of a new 
time scale, and added together, the result is a visual 
picture of the autocorrelation function. 

(9) The process of autocorrelation of a waveform is 
equivalent to passage of the waveform through its matched 
eLlter. 

Figure 2 shows the autocorrelation function of some 


typical waveforms. 
Be CROSSCORRELATION FUNCTIONe® [{Ref. 1,2,3]} 


Autocorrelation can be easily visualized because it can 
be related to the power density spectrun of Fourier 
Analysis. But crosscorrelation does not have a similar 
analogy. 

Crosscorrelation is concerned with the relationship 
between two different signals generated by some common 
process. In general, the crosscorrelation function of two 


waveforms, ese) and ay eel is a graph of the similarity 
between re) and the delayed Bal x Ge) as a function 


of the delay "u" between them. 


Mathematically, crosscorrelation is expressed as: 
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T 
R = ]i V7 Zi .14 
7 A im 1/ _t) x 4(t) x, (ttu) du (251) 


Figure 3 shows a device, Similar to that shown in Figure 
1, designed to obtain the crosscorrelation function fron 
its mathematical definition. In this case, one of the two 
Signals is delayed, multiplied by the other signal and the 
result averaged over a long enough period of time. The 
result is the crosscorrelation function provided that u has 
also been varied over the process. This is also a 


time-delayed crosscorrelator. 
1. Properties of the Crosscorrelation Function. 

The autocorrelation and the crosscorrelation 
functions have different properties. The properties of the 
mmoscecorrelet:on functiam are: 

(1) The crosscorrelation function 1S not an even 


function of u. In other words: 


R (u) #8 = (-u) (2.15) 
Xy xy 


However 


R ({-u) = R- (u) (2216) 
Xy yx 


and this relationship is useful for obtaining R (u) for 
XY 


negative delays. 

(2) The coherence of a signal generated by a 
physical process approaches zero very fast as u approaches 
infinity due to the presence of noise and the uncertainty 


pemniciple. If, added to this, cel and a come fron 


two different, unrelated processes, then: 


18 





‘UBIZBTT YOOTA Joystetioossodg %€ sinzTy : 








oT te "nN sA (N)Y 











SUTBBLOAY fetdsta 


Tattd Fan 


O (9)¢x 


(a) bx 
EE 


19 





R (uu) = OQ (2:00 oa 
xy 


(3) In the frequency domain there exists a Fourier 


transform pair: 


R <=> W -18 
(u) nie (2.18) 


but in this case, oe 1S sometimes called the spectrum of 
Z 
crosscorrelation of the time functions xaOe) and t 5 09) 


Crosscorreiation analysis provides a powerful analytical 


cr 


ool. The ability to measure the degree to which the signals 


cha 


ct 


ise from a2 common physical vhenomenon resemble each 


ry 
ps 


a 
other as a function of the delay time between them, can 
provide a much deeper insight into the phenomenon being 
studied than a separate analysis of the properties of either 


Signal alone. [Ref. 1] 


C. APPLICATIONS OF CORRELATION ANALYSISe 


1. Detection. (Ref. 1,7] 

Correlation analysis plays its most important role 
when noise is to be filtered from very low level signals. In 
this case, noise is any undesired disturbances that mask the 
emaucd transmitted, and in practice, 1t 1S a signal with 
random amplitude variations. In the case of wide band noise, 
the instantaneous value of the signal is nearly independent 
of the value at any other instant which means that the 


coherence time of the process 1s very short. 


20 








Signal detection, such as radar detection, is 
basically the solution of a very important problem common to 
all echo-ranging systems: a signal of known waveform is 
transmitted into a medium and is received, unchanged in fora 
but inmersed in noise. 

When two functions consisting of a signal, S(t), and 


noise, N(t), such as: 


Se c 2 
ee aa ae aU ( ) 


t 3 (Ge N (t 2.20 
Ne) aioe) ee ee ( ) 


are crosscorrelated, the result is: 


-~T 
PG = dia i/27 \" fs (f)4N (t) ILS (t+u) +N (t+u) Jat 
(uj 44 -7.) ( 1 ( 1 ( i 5 5 } 


yA A Tc 

(2.21) 

=R Wu) + uj) + & u) + R u 

ss 'Y sy ys @ yoy 

1 2 eZ 12 1 
(2-22) 
which will approach zero as u approaches infinity. 
Li, xX (ee= eels then the crosscorrelation is: 
1 
R (uy) = R (u) + R- (uy) + R_ (up) + R_ (u) (2225) 
OK SS SN NS NN 


bie in thas case, 
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i 
© 


R = R 
a) (u) 


(2.24) 
NS 


because the signal S(t) and the noise N(t) are uncorrelated. 


ou 1s the autocorrelation function of the 


Original signal and it is a function of u other than zero 


even for large values of u. R (u) is the autocorrelation 
NN 


function of noise. If the transmitted signal is periodic 


then Bat") becomes very small for large values of u_ while 


R (vu) does not. so R (u) is almost equal to R- (u) and 
a0 mx 55 


this result can be used for Signal detection. 
In the case of radar detection, the transmitted and 


the received signals are crosscorrelated. Let 


x, (t) = S(t) + N(t) (2525) 
be the received signal and 
KOC a= eS tt (2226) 
a ) » ' ) 


be the transmitted signal. The frequency of the transmitted 
signal is known so that an internal generator in the 


receiver can be used to generate aCe) as a second input to 


the correlator. When these two Signais are crosscorrelated, 


poe result is: 


— R 22 
yp) Bg TRS (2.27) 


Ze 





Since N(t) and ae are not correlated, then 


R - . 
vs 0 (2225) 


eventually. So 


R =r ; 
Aes (u) (2.29) 


R = 0 (2.30) 


so (u) = 0. This means that if the crosscorrelation 


K 
12 
between transmitted and received signal is not identically 

zero, then a useful signal exists a the received 


transmission. 


2. Direction Finding. [Ref. 3] 

Direction finding is another application or 
crosscorrelation. In this case, the direction of a source of 
arbitrary waveforms (acoustic, electromagnetic, seismic, 
etc), can be determined by crosscorrelating the responses at 
two receivers that may be located far from the source. The 
procedure is represented in figure 4. The waveforms can be 
completely arbitray, and can even be noise. In the figure, 


the wave front gets to receiver 2 first and it will get to 


receiver 1 t seconds later where t = (D/c)cos@ and c is the 
velocity oz propagation. The crosscorrelation has its 
maxinun at a delay time equal to the difference in 
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W Source 


Wave front 


Receiver i , Q 


J = D | aN 


Receiver 2 





X= Deosf} 


ct 
Nt 


(D/c)cosf t Time difference between 
receptions at recelvers. 


Velocity of propagation. 


© 
tH 


Figure 4: Direction Finding. 
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propagation time from the source to each of the two 
receivers. Since D is a constant, then the direction to the 
source, which is the angle 98 measured between the 
propagation direction and the base line between receivers, 
is a function of the delay t (time difference). SO, iia 


source is at o and another is 
at es then the crosscorrelation of the signals at the 


receivers from both sources will peak at two different 
times. 

This technigue can be used to find the direction of 
any disturbance which radiates a plane wave such as cosmic 
noise sources, earthquakes, submarines, and radar or sonar 
targets. 

3. esting Control Systen ~Line. ({kef. 3] 

Tt is often necessarv to determine the transfer 
function or the impulse response of a control system or 
plant which may be in continuous use. [In this case, the 
adjustment should be done without taking the control system 
out of work. Crosscorrelation is idealy used for these 
cases, and the adjustment can be performed on-line. 

Figure 5 representS a typical control syste. 
Low-level wideband noise is introduced into the servo loop 
and it is crosscorrelated with the output. The control 
system acts as a filter and the crosscorrelation represents 
the impulse response of that system. Figure 6 represents a 


typical control system response. 
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Figure 6: Typical Control System response. 





tite GENERATION OF THE CORRELATION FUNCTIONS® 

The correlation Pune e Lens. autocorrelation and 
erosscorrelation, can be generated by three different 
methods. The first two methods (time-domain definition and 
indirect frequency-domain) will be briefly explained in this 
section. The third method (method without pure tine delay) 
Will be explained in the follwing sections. This last 
method can generate the correlation functions in a fairly 
easy and inexpensive way and it is the main concern of this 


thesils. 


Ae HETHOD I: TIMS-DOMAIN DEFINITIONe® [Ref. 1] 


This is an expensive but widely used method. Here, the 
correlation functions are generated by a simple application 
of the mathematical definition of autocorrelation and 
crosscorrelaticn. 

Figure 1 shows a simple device to obtain the 
autocorrelation function and it consists ina step by step 
solution of eguation (2.01). Figure 3 shows the same 
application but for the crosscorrelation function. 

The method consists in delaying one of the signals ‘u" 
units of time, then both signals under consideration are 
multiplied together and the product is fed into a low-pass 
filter. The filter output is one point of the correlation 
function. The complete correlation function is qenerated 
when the delay between the two Signals is varied. Note that 
if the two functions being correlated are equal except for 
the delay "u", the result is the autocorrelation function. 
Otherwise the result is the crosscorrelation function. 

There are some disadvantages to this method, mainly: 

Mineit the two signals are fluctuating, and if the delay 
fu*® is too long, a distortion is introduced in the systen. 


ey) Since this method is a discrete process, the total 
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number otf points which generate the complete correlation 
function depends upon the number of individual delay devices 


used. 


pe METHOD IIT: INDIRECT OR PREQUENCY-DOMAIN METHOD®e 


One of the properties of the correlation functions is 
the fact that they are related to the power spectrum through 
the Fourier Transforms, the autocorrelation being the 
inverse transform of the auto- spectrum and the 
cross-correlation the inverse of the cross-Spectrum. Thus, 
the indirect method consists in transforming to and from the 
frequency domain so that the correlation in the time domain 
is eguivalent to a complex conjugate multiplication of the 
Signals spectra in the frequency domain. 

This method is expensive and its application is based on 
the Fast Fourier Transform algorithm which is a descrete 
(digital) process and 1s mainty used for special purvose 


machanes. [Ref. 1,8] 


Ce WETHOD III: METHOD WITHOUT PURE TIME DELAYe 


The third method used to generate the correlation 
functions is based on the fact that either autocorrelation 
or crosscorrelation can be represented as a Series expansion 
of orthogonal terns. The problem is now how to obtain the 
coefficients of the series expansion and how to use them in 
order to get practical results in a fairly inexpensive and 
Simple way. This method will be fully discussed in the 


foliowing sections. 
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i¥Ve CORRELATION ANALYSIS WITHOUT TIME DELAY 


SE ey ee A ee te —— ae See SS eee ee oe > ay eee ~oee <a i ee a eae ee es ee 


Ae THEORY* [Ref. 2] 


The autocorrelation function can be expanded in a series 


of orthogonal functions as follows: 
CO 
= J 
R(u} = ase (yet pu | OSuc@ (4.01) 
n=0 ona 


where 


a >: coefficients of the series expansion. 
n 


p(u) : a weight function to be defined later. 


8 (u): polynomials that form an orthonormal set with 


respect to p(u) in the range OSu<oO . 


The coefficients a are given by: 
n 


OO 1-g 
a = Jame {u){p(u) ] du (4.02) 
n SB n 


> a ap soy ee ee Se ee 


Figure 7 shows a system that can be used to generate 


the coefficients a. 
n 


If it is supposed that the linear network has an 


impulse response h (t), then the filter output is: 
n 
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co 
xX (n;t) = \x (t-v)h (v) dv 
1 oO 1 nh 


and the multiplier output is: 


CO 
ee) zs) = : \x, (t-v) = (t) ps 


The Output of the low-pass filter is the 


average vaiue of equation (4.04) So: 


cO 


s (n;t) a (t) = - \z {t-v) i (t) Ne (v) dv 


aD 
= een h (v) dv 
O n 


By comparing equations (4.02) and (OG) 


CO 


I-g 
a = Jecne Gi) (2 (eed 
n O n 


and 


CO 


Sans ©) Xer(t) = \rons (v) dv 
1 Z % n 


it can be seen that 
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(4.03) 


(4.04) 


mean Or 


(4.05) 


(4.06) 


nanely: 


ERO 2 


(4.063 








a = xo ne ME) (4.07) 


if h (t) is chosen to be 
n 


1- 
h (t) = 8 (t){p(t)] ° (4.08) 
a n 


Equation (4.08) is a definition. Any filter whose 
impulse response satisfies equation (4.08) is defined as an 
orthogonal filter. 

The above development shows that at least 


mathematicaliy, the coefficients a of the serieS expansion 
n 


can be obtained. Actually, it will be shown next that the 


orthogonal filters can be realized and the coefficients a 
n 


can be obtained from 2 practical system which uses a finite 


humber otf these filters. 
2. Symthesis Problen. 


ee ee PP — = SS eee ee 


Equation (4.01) 


CO 
R(u) = ys a 8 (u){ p(u) a Osu< (4.01) 
=Q nn 
Game be written as: 
2g-1 = lag 
R(a) = [pt] & a6 (u[p(r)] (4.09) 
n=O ono 
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OL 
@ 


2g-1 
R(t) = [p(t) ] : Dak Ct) 0<t<co (4.10) 


where 


1~ 
h (t) = © (t)[p(t)] - (4.08) 
Nn n 


for an orthogonal filter. 
In a practical system, only a finite number of 
filters is used. If N filters are used then equation (4.09) 


can be approximated as; 


i 


ZO = : 
R(t) = [p(t) ] ah (t) (4.11) 
N n=O non 


If g=1/2, which is called the "symmetry" case, then 


“5 
it lz 


R (t) = ah (t) (4.12) 
N Onn 
ama if g = O then: 
N 
a 
R(t) =(p(t)] & ah (t) (4.13) 
N n=0 non 


Figure 8 shows a basic synthesis of the systen. The 
approximate correlation function would ke seen as a 


transient following the application of the impulse. 
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Figure 8: Correlator System without Time Delay 
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A COMVenient system gesilts if in fact the jpulsing 
can be carried out repetitively at intervals that are long 
compared with the filter time constants, and yet 


sufficiently frequently to give a steady trace on a C&T 
screen. [Ref. 1] 


Be TWO REALIZABLE SYSTEMSe [ Ref. 1,2] 

By properly choosing the weight function p(t) and the 
parameter "g" two systems can be realized so that the 
impuise responses of the filters satisfy the definition of 


orthogonality as expressed by equation (4.08). This is so 


because h (t) and 8 (t) depend only on p(t) and "gi. 
n n 


For both systems, the weight function was chosen as; 
=at 
p(t) = ae (4.14) 


fem oystem Lf: q-0. 


When g=0, eguation (4.10) becomes: 
CO 
- 7 3 
R(t) = [p(t)] ah (t) (Gas) 
n=0 nn 
where h (t) reduces to: 
n 


hentia = OG ie) e (4.16) 
n n 


The Laplace transfer function of an kC low-pass 


Peecer 1s: 
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a 
H =—————— = 
_ (3) ae L{ h(t) ] (ae 


where 


a= Vt, = 1/RC (4.18) 


Now 


1 
; HH = ; 
L male Hoe (4.19) 


tt 
mv 
© 


For a simple CR high-pass filter, the transfer 


hunetion is: 





S 
BEG) ee (4.20) 


Now if one low-pass cell and n high-pass cells are 
cascaded together but isolated by a buffer amplifier, the 
resultant circuit will kave a transter function in the 


S-domain expressed as: 


H Fe (4.21) 
N 


Figure 9 shows such a circuit. The inverse Laplace 


transform of equation (4.21) is: 


37 





QD & & Am 


°9=S3 JOJ LEAT TI TeuoVoyuAdo 36 sinsTy 


Momoheim 





38 





teed 
~—-te (4.22) 





h (t) =a 
n 


n 
diz n 


muech@ean be’ written as 


n 
acts kK 
ae (t) = ae dal nl /(n-k)| ii (-at) yx! x! 3 (4.23) 
= ac 
= ae feat) oD (eyelid t) (4.24) 
n n 


where the product is a set of Laguerre functions and lL (at) 
n 


1s defined as a set of Laguerre polynomials. 
It can be seen that equation (4.24) is of the same 


form of eguation (4.08) where L = 8 and so the circuit of 
n n 


figure 9 is by definition an orthogonal filter and can be 


Mseaq tO realize the correlator. 


2. system Ii: gil/2. (Symmetry). 


The realization of this filter is presented briefly 
here for completness Since the experinental work was carried 
out with filters of the first kind (g=0). 

For g=1/2 equation (4.10) reduces to: 


CO 
R(t) = » ah (t) (4.25) 
n=O on 
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and h (t) reduces to: 
n 


1/2 
h (t) = @ (p(t) ] 
n n 


Figure 10 shows the circuit diagram of a 


that can also satisfy equation (4.08). In this case: 


where 
af/2 = 1/RC 
.and its inverse Laplace transform 1s: 


-af2t 
ae 
h(t) = ———— L a 


(4.26) 


filter 


(4.27) 


(4.28) 


4.29) 


which is identical to equation (4.24) except for a constant 


fae tol . 
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The practical work was carried out in two babts. “Firrse- 
the correlator system shown on figure 11 was simulated in 
the computer using the DSL program, and second, experiments 
were conducted on an actual correlator built in the 
laboratory. 

In both cases, the results of the autocorrelation of a 
Sine wave of 155 hz were compared. The practical results 
obtained were very much as predicted by the simulation. 

A detailed explanation of both the simulation and the 


experiment will be shown in the following pages. 


A. NATURE OF THE PROBLEe 


Figure 11 represents the block diagram of the proposed 
weerrelatcer. Powe tne seeasce, “econonv. availability of parts, 
the accuracy of the results and the area in which the 
correlator could be used as a practical device, were the 
guidelines for the design and construction of the correlator 
shown. 

Minimum cost was one of the main objectives of the 
design work. For this reason, the orthogonal filters were 
designed without the buffer amplifiers but this requirement 
limited the range of values of the resistors and capacitors 
necessary to realize the linear filters. Also, the values 
of the resistances and capacitances used were fixed which 
presented the problem of restricting the correlator to 
operate in avery narrow band of frequencies. Fortunately, 
Since the filters are linear, the possibility exists of 
increasing the frequency of operation by varying the time 
Comseant t = RC. 

The low-pass filters, shown after the multipliers, were 
designed as simple RC circuits. 


All of the above considerations mean that several 
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Figure i1: Block Diagrem of practical Correlator. 





non-linearities were introduced in the system. Nevertheless 
the final results were consistently accurate for practical 
purposes. 

The correlator was also restricted to operate at low 
frequencies which together With its narrow-band operation 


made it useful in underwater detection. 


The next problem that had to be worked out was to 
determine the frequency of operation of the train of puises 
necessary to produce the transient needed to generate the 


Seerkricients a. ALtCor some considerations on the 
n 


limitations of the correlator it was decided to use the same 
frequency of the sine wave. The pulses then were in phase 


With the positive peaks of the wave. 


Be LINEAR-FILTER DESIGNe 
Tt can be seen in figure 11 that the total number of 
e 


n 
filters used to generate the coerficients a was equal to 6 


n 

BOEwtne rw OLaGeCrS nh = 0,1,.-..5. Mathematically, an infinite 
hunber of linear filters are required to generate the 
coefficients, but the results obtained by using only six 
filters were very much within the theoretical considerations 
and computer predictions. | 

Higher order fiiters were not used because of the 
availability of components and the characteristics of the 


available pulse generator and oscilloscope. 


Buffer amplifiers have ideally infinite input 
impedance, zero output impedance and unity gain. Figure 9 
shows ah orthogonal filter which uses buffer amplifiers to 


Psolate tne cells. 


dy dy 
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Figure 12 shows a method which can also effectively 
isolate the cells and at the same time eliminates the buffer 
amplifiers. The method consists in introducing a factor of 
isolation k which multiplies the values of the resistors and 
Capacitors in such a way as to maintain the time constant 
the same throughout the cells. [Ref. 1] 

The maximum value of k that can be used depends upon 
the number of celis and the available values of resistors 
and capacitors. It was found that for this waork, a factor 
of isolation of over 10 would require values much too big 
for the resistors and too small for the capacitors. For k 
less than 5, the approximation of the ideal impulse response 
was not accurate enough. A value of 5 was finally used and 
the results are presented in figures 19,20, and 21. 

The number of celis of the filters depends upon the 
degree of the Laguerre polynomial. The low-pass R&C cel 
represents the Laguerre Polynomial of degree 0 (n=0). Th 
degree is increased by cascading high-pass RC cells. The 
mhaxinua number of cells that could be used without 
introducing appreciable distortion was found to be 6 {(n=5): 


one low-pass cell and five high-pass cells. 
2. Details of circuitry 


A simplified schematic diagram of the orthogonal 


filter for n=5 is shown in Figure 13. Theoretically: 


and the following values were chosen: 
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R = 1000hn C = 10 MP 
0 0 


a= 1/RC = 1000 k 


ul 
wn 


Table I shows the actual values used together with 
the theoretical values. 

The filters required for the lower orders can be 
obtained by simply eliminating one high-pass cell for each 


order. Figures 19, 20 and 21 shows the actual results. 


TABLE [Le 
Theoretical Values Actual Values. 
R C R K C K a 

R Cc 

R wtk ¢ 10MF - iK 10MFP 1000.00 
0 0 

ae ao eS MF 5.1K soe Ze Sie) 9890.39 
225K > ~4YOMF Ze UKs eHIME 4.6 1073.17 
sa 2 cK oo -OSME 122. 51K 520 “S0S2NEF 8.0 974.83 


R Ogio oF MOT oul Wes) tl 85 10 Oe My ue 7 oes 


ie S02. OK c. -0032MF SI 2K 0h. Cee UUs otha. | 931.46 


Ce MULTIPLIERS®e 


Integrated circuit multipliers were used to realize the 
correlator. In this case, the Intersil] 8013A was used which 
is a four quadrant analog multiplier and whose output is 


Proportional to the algebraic product of two input signals. 
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Fisure 14: Block Diagram of Multiplier. 
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Figure 15: Connection Diagram of Multiplier. 
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The multipliers can operate with standard + 15 volts 
Supplies, and have a maximum dissipation power of 500 mu. 

The block diagram of the multipliers is shown in figure 
14 while figure 15 shows the connection diagram together 
With the trimming circuit necessary for gain accuracy, 


offset voltage and feedthrough performance. 


De LOW-PASS FILTER DESIGNe [{Ref. 4} 


The low-pass filters were designed as _ siaple ie 
low-pass cells with a cut-off frequency of about 600 dz. 
More accurate filters can be designed, such as second order 
Butterworth filters, but in this case, the simple R-C 
low-pass filters were accurate enough for practical 
purposes. 


The designed cf the filters was catried out as follows: 


oo 2 (5.01) 
2RC 





If R is chosen to be 510 Ohms and C equal to 0.47 NF, 





then the cut-off frequency would be: 


[7 nani een 
-~6 
2 (510) (0.47x10 ) (5.02) 


For the computer Simulation, a real-pole transfer 


function is required. This transfer function is given by: 





1 

H S = =e 2 

ape ) (5.0 3) 
Ps+1 RCs+1 


-6 -4 -4 
where P= RC = 510(0.47x10 )= 2.397x10 = 2.4x10 


and P iS 2 parameter required for the real-pole transfer 
m 


function computer block. 


Ee THE ADDER ANPLIFIERe [Ref. 4] 

The output of the adder or summing amplifier of figure 
17 is a linear combination of the input signals. This 
arrangememt was used to construct the adder required in the 
last stage of the correlator. 

As will be explained in the next section, the output of 
the amplifier is not the true correlation function but the 
product of the correlation with the weight function used to 
realize the linear filters. Nevertheless, the shape of the 
output can very well be used for signal identification since 
it is a characteristic of the signal being correlated. 


Figure 18 is the circuit diagram of the correlator. 
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Figure 17: Diagram of Adder Amvlifier. 
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Figure 18: Correlator Circuit Diagram. 
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Figure 19: Imoulse Hesvonse of Linear Filters, 
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Fizure 20: Imoulse HKesvonse of Linear rilters. 
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Fizure 21: Imvulse Kesvonse of Linear Filters. 


n=5 
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tiszure 22: Practical Correlator Output. 
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VIe COMPUTER SIMULATIONS 

Figure 22 shows the diagram used to simulate the 
correlator on DSL (Digital Simulation Language) together 
With the name of the variables used for the program which 

appears in appendix A. 

| Since the linear filters used to build the correlator 
were those for which g=0, the actual output of the 
correlator 1S not the true correlation function but the 
correlation function multiplied by the inverse of the weight 
Pinckason 9p Ct). In the case of the autocorrelation of the 
Sine wave, the output of the correlator which iS shown in 


Figure 24, 1s: 


R' (u) = ae R (u) (6.01) 


In order to get the true autocorrelation function of the 


Sine wave shown in Figure 23 R' {u) has to be multiplied by 
KX 


at 
the function e /a shown on figure 25 This function is an 


exponential which increases very rapidly and overflows the 
computer results for times beyond those shown on the figure. 
Scaling can be used to compensate the overflow and increase 
the time of observation of the autocorrelation function. 


Thus: 





R (u) = R' (u) (6.02) 
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Pigures 28 to 33 show the response of the linear filters 
for a train of impulses. They can be compared to Figures 19 
to 21 which show the impulse resvonse of the filters. The 
pictures were obtained from the experimental work. 

The response of the filters are of the same form as the 
Laguerre Functions shown on Figure 27 and which are 


necessary to generate the coefficients a of the series 
n 


expansion. 

Figure 22 shows the output of the experimental 
correlator. This figure can be compared to Figure 25 which 
represents the output of the simulation. It can be seen 
that the experimental results were in agreement with the 


computer predictions. 
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VII® CONCLUSIONS® 


The orthogonal filters were designed in such a way as to 
eliminate the buffer amplifiers which were required to 
isolate the cells. This requirement, together with the fact 
that the values of the resistances and capacitances’ used 
were fixed, restricted the correlator to operate in a narrow 
band of frequencies. 

Some work can be done in this area in order to 
investigate the possibility of increasing the range of 
frequencies of operation of the correlator by means of 
variable resistances and capacitances which will vary the 
time constant t=RC. This is entirely possible because the 
filters are linear. 

Also, further investigation may be reguired for the use 
of some other type of multipliers in such a way that they 
Can be grouped aS one unit instead of uSing 12 separated 
units as designed. . 

The actual output of the correlator is not the true 
COEFEGlatwoneefunMccion but 2S an exponentially decayed 
correlation as shown in Section VI. So, it would he worth 
investigating on the possibility of obtaining the true 
correlation function by the generation of an exponential of 

fact 


the type e fa in order to compensate the exponentially 


decayed correlations obtained in this study. 
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